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ABSTRACT 

■ — A  closed-form  solution  for  the  Timoshenko-type  shear  correction 
coefficient  (KiJ  governing  the  deflection  of  bimodular  composite  laminates 
in  cylindrical  bending  is  presented.  The  bending-stress  distribution  for 
a  laminate  constructed  of  bimodular  materials  (which  have  different  elastic 
moduli  in  tension  and  compression)  is  used  in  the  two-dimensional  equili¬ 
brium  equation  to  obtain  the  transverse  shear-stress  distribution.  This 
shear-stress  distribution  is  used  to  obtain  expressions  for  the  shear  cor¬ 
rection  coefficient  (based  on  equivalent  shear  strain  energy)  and  the  maxi¬ 
mum  dimensionless  transverse  shear  stress  (i  )  .  Finally,  the  effects  of 

xr  max 

the  elastic-constant  ratios  on  the  neutral -surface  location  and  shear  cor¬ 
rection  coefficient  for  laminates  consisting  of  either  unimodular  or  bimodular 
materials  are  studied.^”  / 

(Tau  /  v!  ^  x 


INTRODUCTION 


Materials  which  have  different  moduli  in  tension  and  in  compression 
are  called  bimodular  materials.  Rock,  concrete,  cord-rubber,  paperboard, 
and  certain  biological  tissues  are  examples  of  such  materials.  Even  aramid- 
fiber,  polymer-matrix  composites  exhibit  some  bimodularity.  The  analysis  of 
laminated  bimodular  material  is  more  complicated  than  unimodular  material 
(ordinary  material)  due  to  the  dependency  of  the  material  stiffness  on  the 
material  properties,  which  indeed  depend  on  the  state  of  stress  ( i . e . , 
tensile  or  compressive)  in  the  laminate. 

Although  transverse  shear  deformations  have  been  considered  in  the  analyses 
of  bimodular  laminates  [1-4]  in  recent  years,  there  has  been  no  effort  to  in¬ 
clude  the  effect  of  bimodularity  of  the  material  on  the  Timoshenko- type  shear 
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correction  coefficient,  K  .  In  other  words,  it  has  been  tacitly  assumed  that 
2 

the  value  of  K  is  the  same  as  that  of  an  ordinary-material  laminate.  The  use 
of  such  a  shear  correction  factor  in  predicting  the  deflection  of  shear  deform¬ 
able,  ordinary-material  laminates  in  cylindrical  bending  is  well  established 
[5,6].  It  is  noted  that  Whitney  [5]  extended  Chow’s  symmetric- laminate  work 
[7]  to  arbitrary  laminates.  Also,  it  can  be  shown  that  the  resulting  expres¬ 
sions  in  [5,6]  are  algebraically  equivalent. 

In  this  paper,  a  straight-forward  approach  analogous  to  that  used  in 
elementary  shear  theory  for  single-layer  ordinary  materials  is  employed. 

It  may  be  considered  to  be  a  generalization  of  the  work  of  Bert  [6]  from 
ordinary-material  laminates  to  bimodular-material  laminates. 

THEORY  AND  FORMULATION 

Consider  a  rectangular-cross-section  laminated  beam  of  thickness  h.  The 
origin  of  the  Cartesian  coordinate  system  is  located  on  the  mid-surface  of 
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the  beam  with  the  x  axis  being  along  the  length  of  the  beam  and  the  z  axis 
being  measured  positive  downward.  The  same  displacement  field  as  used  in 
shear  deformable  beam  theory  is  implemented  here: 

u(x,z)  «  u° ( x )  +  z v ( x ) 

o  (1 

w(x,z)  =  w°(x) 

The  axial  normal  strain  and  transverse  shear  strain  are  given  by 

S  =  U!x  *  Zi,X  '•  YXZ  ’  •  +  “!x  <2 


where  a  comma  denotes  the  derivative  with  respect  to  the  quantity  follow¬ 
ing  it.  The  longitudinal  bending  stress  at  any  distance  z  from  the  mid¬ 


plane  of  the  laminate  is 


=  E(k)s 


k=c,t 


t  =  G(k)- 


The  longitudinal  stress  resultant  and  stress  couple  are  defined  as 

h/2 

(N,M)  =  |  (l,z)oxdz 
-h/2 

With  the  aid  of  Equations  (2),  (3),  and  (5),  one  can  derive  the 
following  laminate  constitutive  relation 


’a  8 lk, 
_8  °Jk> 


where 


I  1/  «- 

(A,B,D)  =  f  E(k)(l,z,z2)  dz 


k=c,t 


In  the  absence  of  body  force,  the  two-dimensional  static  equilibrium 
equation  for  forces  in  the  direction  along  the  length  of  the  beam  is 


O  +  T  =0 

x,x  xz.z 


(7) 


Integrating  Equation  (7)  with  respect  to  z  and  using  Equations  (3)  and  (4), 
one  has 


-  -  [G(k)(A0-  B2)!1  C ( Da  -  Bb ) N  +  (Ab-Ba)M  J 

aZ  |  a  y  A 


(8) 


where  a  and  b  are  "partial  stiffnesses"  for  stretching  and  bending¬ 
stretching  coupling  defined  by  the  following  expressions  for  the  convex 
bending  case  (top  layers  in  compression  and  bottom  layers  in  tension) 


(l,z)Ec  dz 


(a,b)  = 


-h/2 


z„  >  z  >  -h/2 
n 


(9) 


(1 ,z)Ec  dz  +  (1 ,z)E  dz  h/2  >  z  >  z. 


-h/2 

Multiply  Equation  (7)  by  unity  and  z,  respectively,  and  then 
integrating  through  the  thickness,  one  obtains  the  following  equilibrium 
equations,  which  coincide  with  those  of  elementary  beam  theory. 


M  ■  Q  ,  N  -  0 

•  *  t  A 


where 


h/2 

Q  -  yxz  f  G<kl  dz 
-h/2 

Therefore,  Equation  (8)  simplifies  as  follows: 


(10) 


(ID 


J 

U 

M 


8 


L; 

V] 

’.i 

fa 
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(12) 


Yxz  *  -  Q(Ab-Ba)LG(k)(AD-B2]-1 

Extending  the  definition  of  the  shear  correction  coefficient  as  used  in 
[6],  which  is  based  on  the  shear  strain  energy,  to  the  case  of  bimodular 
material,  one  has 


K2  - 


(AD-  BZ) 


h/2 


.(k) 


h/2 


(13) 


dz][ 


-h/2 


(l/S(k))(Ab-  Ba)2dz] 


-h/2 


If  G  is  constant  throughout  the  thickness  of  the  beam.  Equation  (13)  can 
be  simplified  to 


K2 


hf 


2  2 

(AO  -  B  j 
h/2 

(Ab-Ba)2 


dz 


(14) 


-h/2 


By  changing  the  coordinate  from  z  measured  from  the  midplane  to  z'  measured 
from  the  neutral-surface  position,  one  can  simplify  Equation  (14)  as  follows: 


K2 


(15) 


To  evaluate  the  stiffnesses  (A, 8,  and  D),  one  needs  to  know  about  the  loca¬ 
tion  of  the  neutral -surface  position.  A  closed-form  solution  can  be  obtained 
only  when  the  neutral -surface  position  (zn)  and  thus,  the  stiffnesses  do  not 
depend  on  x.  The  same  criterion  to  define  zp  as  used  in  [4]  is  applied  here 

cx  *  0  (16) 
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Using  Equations  (2)  and  (5)  and  assuming  that 


one  gets 


N  =  0 


zn  =  B/A  =  constant 


(17) 

(18) 


Equation  (18)  can  be  solved  explicitly  for  z  . 

n 

It  is  convenient  to  introduce  the  dimensionless  transverse  shear 
stress  defined  as  follows: 


xz 


xz 

Q/h 


(19) 


Combining  Equations  (4)  and  (12)  and  measuring  z1  from  the  neutral- 
surface  location,  one  obtains 


*xz  =  bh/0 


(20) 


Detailed  derivations  for  K  and  (f,J  are  given  in  the  Appendix. 

xz  max 


NUMERICAL  RESULTS 

Three  cases  of  laminates  constructed  of  bimodular  materials  are  con¬ 
sidered  here:  single-layer,  two-layer,  and  three-layer.  Numerical  results 
obtained  from  closed-form  solutions  for  the  shear  correction  coefficient 
and  the  maximum  dimensionless  shear  stress  for  different  bimodular-material 

parameters  are  presented.  Also,  the  effects  of  elastic-constant  ratios  on 

2 

neutral -surface  position  and  K  for  both  unimodular  and  bimodular  laminates 
(up  to  three  layers)  are  studied.  The  laminates  are  assumed  to  be  made  of 

equal -thickness  layers  and  the  effect  of  the  sign  of  the  longitudinal  normal 

* 

Note  that  this  assumption  is  made  only  for  the  bimodular  case  (not 
unimodular) . 


c  t 

strain  on  shear  moduli  is  ignored  (that  is*  G  / G  =1). 

In  the  case  of  single-layer  bimodular  material  (see  the  Appendix),  it 

is  interesting  to  note  that  both  the  transverse  shear  correction  coefficient 

and  the  maximum  transverse  shear  stress  are  unaffected  by  the  bimodular 
c  t 

ratio  E  /E  .  However,  this  does  not  imply  that  the  distribution  of  shear 

stress  is  the  same  for  unimodular  and  bimodular  materials  (see  Fig.  A-l). 

The  effect  of  bimodular  ratio  on  the  location  of  the  neutral  surface  is 

c  t 

studied  in  Table  1.  For  values  of  E  <  E  (e.g.,  aramid-cord/rubber) ,  h 

falls  within  the  lower  half  of  the  thickness;  whereas  for  materials  such  as 
c  t 

rock  (E  >  E  ),  h^  is  in  the  upper  half  of  the  beam  (for  the  convex  bending 
case) . 

In  bending  of  a  two-layer  bimodular  laminate  since  one  of  the  two  layers 

is  always  in  compression  (or  tension),  only  three  of  the  four  elastic  moduli 

belonging  to  these  layers  are  pertinent  [see  Fig.  A-2].  This  means  that 

the  location  of  the  neutral  surface  depends  on  only  two  elastic-moduli 

ratios.  The  effects  of  these  two  ratios  on  ht  are  studied  for  a  wide  range 
c  t 

of  E2  /E 2  in  Table  2.  The  shear  correction  coefficient  varies  drastically 

with  the  ratio  of  the  shear  moduli  of  the  two  layers  (G-i^),  whereas  it 

r  t  c  t 

changes  only  a  little  for  a  wide  range  of  E2  f°r  E ]  =  1  *  as  shown 

in  Fig.  1 . 

The  shear  correction  coefficient  for  two-layer  unimodular  laminates 

increases  rapidly  from  50%  to  98%  of  the  classical  value  (5/6)  as  the  elastic- 

modulus  ratio  E^/Eg  varies  .from  0  to  0.25  as  shown  in  Fig.  2.  For  ratios 

2 

greater  than  0.5,  a  value  of  5/6  would  be  a  good  approximation  for  K  .  The 
dimensionless  maximum  shear  stress  increases  from  0  to  12  as  E^/Eg  changes 
from  0  to  infinity  as  tabulated  in  Table  3. 
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Three-layer  bimodular  laminates  which  have  facings  (top  and  bottom 

layers)  made  of  the  same  material  are  considered  here.  For  the  case  in 

which  E-j0/^0  =  1,  the  effect  of  bimodular  ratio  on  h ^./h  for  differ- 

ent  ratios  of  2^  (1/2,  1,  and  2}  are  investigated  in  Table  4.  As 

E^/E^  increases  from  0  to  10,  h^/h  decreases  for  all  three  values  of 

^2  ^2 '*  ^so’  a  ^  vs  E^/^  1S  9''ven  ^9*  3  which  shows  the 

shear  correction  factor  as  a  function  of  the  bimodularity  of  the  material. 

2 

Finally,  the  affect  of  elastic  moduli  ratio  E9/E,  on  K  and  (t  ) 

c  1  xZ  max 

are  studied  (see  Fig.  4  and  Table  5)  for  three-layer  unimodul ar-material 

2 

laminates  with  identical  facings.  For  different  values  of  G-j/G^,  K  increases 
rapidly  for  E-j^  less  than  0.05  but  after  this  point,  K  changes  very  slowly. 
The  maximum  dimensionless  shear  stress  varies  from  0  to  40.5  for  0  <  E^/E^  <  ». 

CONCLUSIONS 

Closed-form  solutions  for  the  Timoshenko-type  shear  correction  coefficient 
and  the  maximum  dimensionless  transverse  shear  stress  are  presented  for  bi¬ 
modular  laminates  undergoing  cylindrical  bending.  These  solutions  depend  on 
bimodular  ratio(s)  for  different  laminations.  However,  it  is  interesting  to 
note  that  both  the  transverse  shear  correction  coefficient  and  the  maximum 

dimensionless  transverse  shear  stress  (x  „)  for  a  single-layer  material 

xz  max 
c  t 

are  unaffected  by  bimodular  ratio,  E  /E  . 

2 

For  two-  and  three-layer  laminates,  both  K  and  (x  _)  depend  upon  the 

xz  max 

bimodular  ratios  and  transverse  shear  moduli  ratios.  According  to  the  results 

2 

presented,  in  some  cases  K  is  less  than  the  classical  value  of  5/6. 
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APPENDIX:  DETAILED  DERIVATIONS 


In  the  following,  detailed  derivations  for  unimodular  and  bimodular 
laminates  are  given.  Three  types  of  lamination  are  considered:  single¬ 
layer,  two-layer,  and  three-layer.  In  the  latter  two  cases,  all  layers  are 
assumed  to  be  of  the  same  thickness.  Also,  in  the  bimodular  cases,  the 
effect  of  bimodularity  on  elastic-shear  moduli  is  ignored  (since  for  most 
of  materials  G  /G  =1).  Further,  in  the  bimodular  analyses,  axial  force  is 
assumed  to  be  zero. 


★ 

1.  Single-Layer  (Bimodular  Material )  ( 

Normal  and  shear  stress  distributions  for  single-layer,  bimodular  mate¬ 
rial  are  shown  in  Fig.  A-l .  Since  the  axial  force  is  taken  to  be  zero. 


N  =  [  ax  dz'  *  <  [j  E  z'dz'  + 

-h_ 


Etz'dz'] 


*  0 


(A-l) 


Using  the  fact  that 


h  =  ht  *  nc 


(A-2) 


one  can  integrate  Equation  (A-l)  to  obtain 


h7h  = 


,  h. /h  = 


/B 


1  +  /B  1  1  +  /s 


(A-3) 


where  3  is  the  bimodular  elastic-moduli  ratio  defined  as 

B  =  Ec/Et 

From  Equations  (6)  and  (9),  one  has 


*This  case  also  might  be  considered  to  be  a  two-layer  laminate  of 
ordinary  materials  but  with  a  certain  special  thickness  ratio 


D  =  (EChc3  +  Etht3)/3 

(A-4) 

b=I(EC(z‘)2-EChc2 

-  h_  <  z  ’  <0 

c 

(A-5) 

2  (  Efc(z ' )2  -  Echc2 

0  <  z'  "  ht 

substitution  of  Equations  (A-4)  and  (A-5)  into  Equation  (15)  and  using 
Equation  (A-3)  leads  to 

K2  =  5/6  ( A-6) 

Also,  using  Equations  (19),  (20),  (A-3),  (A-4),  and  (A-5),  one  has 

(t„J  =  3/2  (A-7) 

max 


Although  the  maximum  shear  stress  remains  unaffected  by  bimodularity  ratio, 
the  shear-stress  distribution  does  depend  on  3.  Combining  Equations  (19), 
(A-3),  (A-4),  and  (A-5)  gives  the  following  shear-stress  distribution 


xz  2 


3  (  /3  -  /3  (/s  +  l)2(z/h)2 

(  Ss  -  —  (/s  +  l)2(z/h)2 
/s 


/s  +  l 


<  z/h  <  0 


(A-8) 


0  <  z/h  < 


/3  +  1 


Two-Layer  Laminates 


(a)  Bimodular  Material 


Each  individual  layer  has  different  elastic  moduli  in  tension  and 
compression.  However,  since  one  layer  is  only  in  either  a  tension  or  a 
compression  state,  only  three  of  the  four  elastic  moduli  for  two-layer 
analysis  come  into  the  picture  (see  Fig.  A-2).  A  similar  procedure  to 
that  used  for  the  single-layer  case  is  used  here.  To  determine  the  bending 
and  the  partial  bending-stretching  coupling  stiffnesses,  one  is  able  to 
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obtain  the  neutral -surface  position  from  the  following  quadratic  equation 


-h  +h/2 


N  =  <[(  C  EjCz'dz‘  +  f  E2cz'dz'  + 


E2tz'dz‘]  =  0 


(A-9) 


-h 


-hc+h/2 


Knowing  the  neutral -surface  location,  one  can  easily  determine  the  follow¬ 
ing  sti ffnesses 

3.i  3 1  u  \  3  i  r  tu  3^ 


D  =  {E1c[(h/2-  hQ)3  +  hc3]  +  E2c(h/2-hc)J  +  E2 \  1/3 


(A- 10 


and 


b  =  A- 


E1C(21 )2 

P  c  h  2 
"  E1  hc 

-hc  <  z 

E2C(z')2 

p  t  .  2 

'  E2  ht 

-h  + h/2  <  z 

E2t(z* )2 

p  t  .  2 

'  E2  ht 

0  <  z 

( A- 1 1 


Substituting  Equations  (A-10)  and  (A-ll)  into  Equation  (15),  one  has 


/2  _  5 


16  gi[(q+s)(h/2-hr)3  +  --th^h^] 


6  h(g-|  +  1 )  [3(g1oi2  -  S2) (h/2  -  hc)5  -  10(gia2h£  -  3h2)  (h/2  -  hc)  : 


+  15(gia2h^  -  hj) (h/2  -  hc)  +  8 ( g-j a2 h^  +  h*)]  (A-12 

where 

a  =  E1c/E2t  ,  3  =  E  2  c/  E  2 1  and  g]  = 

2 

Since  hc/h  and  h t/h  depend  upon  a  and  3,  K  is  a  function  of  g-j ,  i,  and  3 
for  the  case  of  equal -thickness  layers.  The  maximum  dimensionless  shear 
stress  also  is  a  function  of  the  bimodular  ratios  (i,3) 


(  a[h“-(h/2-h  T] 

(-  )  =2 _ h _  1  'hc  z'  '  'hc+h/2 

xz  max  2  g[(h/2-h)3+h3]+B(h/2-h)3+h3  ) 

1  <  V  0  -  z'  ht 

(A-l 3) 

(b)  Unimodular  Material 

A  similar  lamination  as  used  in  part  (a)  is  considered  here.  The  normal 

and  shear  stress  distributions  are  shown  in  Fig.  A-3.  Since  the  flexural- 

extensional  coupling  stiffness  (b)  does  not  vanish,  one  must  use  Equation 

2 

(13)  in  order  to  compute  K  .  Let  the  transverse  shear  moduli  be  the  same 

★  _ 

for  both  layers  (G  /G.  =1)  and  g  =  ^^l"  Then,  the  stiffnesses  can  be 
expressed  as 

A  =  E1h(a-l)/2 

8  *  E-j h2(a-l  )/8  (A-14) 

D  =  E1h3(a+l)/24 
and 

(  C(z  +  h/2)  ,  (z2  -  h2/4)/2]  -h/2  <  z  <  0 

(a,b)  =  E,  {  ?  (A-l 5) 

1  (  [(az  +  h/2)  ,  (aZ-hV4)/2]  0  <  z  <  h/2 


Substitution  of  Equations  (A-14)  and  ( A-l 5 )  into  Equation  (13)  gives 


K2  .?• 


5_  .  _ (a2  +  14a  +  1 ) _ 

^  2[g4  +  1 2g3  +  102g2  +  1 2g  +  1  ] 


( A-l 6 ) 


It  is  interesting  to  note  that  as  given  by  Equation  ( A- 1 6 )  gives  the 
same  result  for  g=c  and  g=l/c,  where  c  is  an  arbitrary  constant.  Combining 


★ 

This  would  be  true  if  the  individual  layers  were  made  of  the  same 
orthotropic  material  but  oriented  at  0°  and  90°,  respectively,  i.e.,  a 
cross-ply  laminated. 
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Equations  (4),  (12),  (19),  (A-14),  and  (A-15)  one  gets 


{'TXz\ 


3  .  8a(  >H) 

^  a2  +  14a  +  1 


(A-17) 


which  shows  that  (t  )  depends  on  elastic-moduli  ratio  (E«/E,),  only. 

xz  max  c  1 


3.  Three-Layer  Laminate 
(a)  Bimodular  Material 

Consider  a  three-layer  laminate  with  the  top  and  bottom  layers 

/•  /»  f  f 

(facings)  made  of  the  same  material  (i.e.,  E^  =  E^  ,  = E3  )  as  shown 

in  Fig.  A-4  in  the  convex  downward  bending  position.  One  assumes  that  the 
neutral-surface  position  is  within  the  middle  layer  (which  happens  for 
most  practical  cases).  The  axial  force  vanishes  if 


-hc+h/3 


EiZ'dz'  + 


E2cz'dz'  + 


-h  +h/3 

V* 


h^-h/3  h 

E2tz,dz  ’  +  [  E^z'dz'  =  0 


ht-h/3 


( A- 18) 


Equation  (A-18)  can  be  used  to  determine  the  neutral -surface  position. 

2 

Having  the  values  of  h^  and  hc  enables  one  to  calculate  K  with  the  same 
procedure  as  used  for  the  two-layer  case 
E  c 

D  =  (a[h^  +  (h/3-hc)3]-  (h/3-hc)3  +  2(ht-h/3)3  +  Y[h^-  (ht  -  h/3) 3]  > 


( A- 19) 


G  k  dz  *  G2h[2g1  + l]/3 


(A-20) 
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i  (G~  )  b^dz  =  -jl -  {{gTa2{i[(h/3-hc)5+hc5]-|hc[(h/3-hc)  +hc]  +  hh^/3} 

ihc  2 

-Cj(h/3-hc)5  +  |  H*(h/3-  hc)3  +  (H2)"(h/3-hc)] 

+^ht"h/3)5  +  T  Hc(ht*  h/3)3  +  (Hc)2(ht"  h/3)] 
+91'1y2{J-ChJ  -  (ht-h/3) 5]  -  |h2  [h2-(ht-h/3) 3]  +  hhJ/3}  }} 

(A-21 ) 

where 

Hc2  =  (h/3"hc)2(3l'1)  “  hc2  (A-22) 

a  =  E1C/E2C  ,  3  =  E2t/E2C  ,  y  =  E3t/E2C,  and  g]  =  G^Gg 


By  substituting  Equations  (A-19),  (A-20),  and  (A-21)  into  Equation  (15),  one 
2 

can  determine  K  .  Also,  using  Equations  (A-19)  and  the  following  expression 
for  bmax  (the  maximum  value  of  partial  coupling  stiffness  b)  in  the  respective 


regions,  one  has  the  maximum  dimensionless  shear  stress  by  Equation  (19). 


11  e  c 
Ll- 1 

L  E  C  -  L  E  C 
>ltl  L2t2 

■  r  C  -  r  C  ,  r  t 

Lr  1  '  L2L2  L3L2 

LylE,1 


-h  <  z'  <  -h  +h/3 

-hc+h/3  <  z'  <  0 

0  <  z1  <  h^-h/3 


hf-h/3  <  z'  <  h. 


[26a  + 1]  g-, 


(l+2g1 )[102a2  +  g1 (120a2  +  20a  +  1 )] 


Absolute  dimensionless  maximum  shear  stress  occurs  at  the  midplane. 
Using  the  same  procedure  as  used  before,  one  obtains 

( Txz)max  =  (tx2)2=0  =  \  '  ^Sa+T  (A' 

and  the  shear  stress  at  the  interfaces  is 

\  _  3  24a 

Txz  7  =  +  h/6  ‘  7  *  26a  +  1  (A‘ 


It  is  interesting  to  note  that 


^Txz^ax 


xz'z  =  +  h/6 


=  1  +  ( 1 /8a) 


The  top  and  bottom  layers  have  the  same  properties. 


SYMBOLS 


a 

A 

b 

B 

c 

D 

E 

9 

G 

h 

Hc 

K2 

L.  1-1.2.3.4 

M 

N 

Q 

S 


u,w 


X 

z 


z' 


z 


n 


partial  excensional  stiffness 
extensional  stiffness 

partial  flexural-extensional  coupling  stiffness 

flexural-extensional  coupling  stiffness 

arbitrary  constant 

flexural  stiffness 

longitudinal  Young's  modulus 

transverse  shear  moduli  ratio 

transverse  shear  moduli 

beam  thickness 

defined  by  Equation  (A-22) 

shear  correction  coefficient 

defined  by  Equation  (A-24) 

moment  resultant 

force  resultant 

transverse  shear  resultant 

thickness-shear  stiffness 

displacements  in  x  and  z  directions 

coordinate  along  the  length 

thickness-direction  coordinate,  measured  positive  downward 
from  the  midplane 

thickness-direction  coordinate,  measured  positive  downward 
from  neutral -surface  position 

neutral -surface  location  for  zero  normal  strain  in  the 
x  direction 
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*>3,v 


bimodular  longitudinal  Young's  moduli  ratios 
normal  strain  in  the  x  direction 
Yxz’Yxz  engineering  shear  strain  in  xz  plane:  actual  and  uniform 
shear-distribution  cases 
i  shear  stress  in  xz  plane 

At 

cr  normal  stress  in  the  x  direction 

bending  slope 

Subscript 

1 ,2,3  layer  number 

(  ).x  3(  )/3x 

Superscript 
o  midplane 

dimensionless 

(c,t)  (compression,  tension) 

(k)  index  (=  c,t) 
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Table  1.  Effect  of  Bimodular  Ratio  on 
Neutral -Surface  Location  for  Single- 
Layer  Bimodular  Material 
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Table  2. 

Effect  of  Bimodular  Ratios  on 

Neutral - 

r> 

Surface 

Location  for  Two-Layer  Bimodular 

Laminates 

« _- 

(h1/h2  = 

1) 

i 
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2  xt2 

ht/h 
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Table  3.  Effect  of  Elastic-Moduli  Ratio  on 
Dimensionless  Maximum  Shear  Stress  for  a  Two- 
Layer  Unimodular-Material  Laminate  (h^/h2  =  l) 
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Table  4.  Effect  of  Bimodular  Ratios  on  Neutral- 
Surface  Location  for  Three-Layer  Bimodular  Lamin¬ 
ates  (h]/h2*h3/h2  =  1,  E^/EgC  =  E3c/E2c  =  1 ) 
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Figure  A-2.  Normal  and  shear  stress  distribution  in 
a  two-layer  bimodular  material  laminate. 
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